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Abstract 



Context. To understand the feedback of black holes on their environment or the acceleration of ultra-high energy cosmic rays in the 
^■"^ present cosmic epoch, a comprehensive inventory of radio galaxies in the local universe is needed. This requires an all-sky catalog of 

■ ^ radio-emitting galaxies, that hitherto has not been available. 

* Aims. We present such an all-sky sample. Our catalog allows one to build volume-limited subsamples containing all low-power radio 
galaxies, similar to the prototypical low-power radio galaxies Cen A or M87, within some hundred Mpc. 

CN Methods. We match radio emission from the NVSS and SUMSS surveys to galaxies of the 2MASS Redshift Survey (2MRS) using 

an image-level algorithm that properly treats the extended structure of radio sources. 

Results. The bright master sample we present contains 575 radio-emitting galaxies with a flux greater than 213 mJy at 1.4 GHz. 
Over 30% of the galaxies in our catalog are not contained in existing large-area extra-galactic radio samples. We compute the optical 
f) and radio luminosity functions and the fraction of radio galaxies as a function of galaxy luminosity. 94% of the radio galaxies within 

• z = 0.03 are of Hubble type E/SO. The local galaxy density in a sphere of 2 Mpc centered on the radio galaxies is 1.7 times higher 
than around non-radio galaxies of the same luminosity and morphology, which is a statistically significant enhancement (> 3cr). 

Oh Conclusions. Our sample presents the deepest all-sky catalog of low-power radio galaxies. The observed enhancement of the galaxy 

' density around radio galaxies suggests a causal relation between external galaxy properties, such as environment or merger history, 

^ and the formation of powerful jets in the present universe. Since the enhancement is observed with respect to galaxies of the same 

i luminosity and Hubble type, it is not primarily driven by black hole mass. Our automated matching procedure is found to select 

^ radio-emitting galaxies with high efficiency (99%) and purity (91%), which is key for future processing of deeper, larger samples. 



> 1. Introduction 

(T^ The final episode in the history of black hole accretion and 
galaxy formation takes place in our cosmic backyard, the lo- 
cal universe. A large sample of nearby radio galaxies may thus 

\^ be considered an important anchor point for theories of black 
hole growth and downsizing (e.g., Alexander & Hickox 2012) 

fSJ and can be used to study feedback from radio jets (Fabian et al. 

^— H 2003; Mathews & Brighenti 2003) on the environments of their 
^ host galaxies. 

• ^ Nearby radio galaxies allow one to study a regime of fainter 
jets and thus, assuming standard jet-disk coupling (Rawlings & 
J-H Saunders 1991; Falcke & Biermann 1995), lower rates of accre- 
^ tion onto the black hole. At these low rates, the mode of accre- 
tion (Narayan & Yi 1995) may switch to a 'jet-dominated mode' 
or 'radio mode' (Falcke, Kording, & MarkoflT 2004), similar to 
the 'low-hard state' of X-ray binaries (Fender, Belloni, & Gallo 
2004), implying that radio surveys are an ideal tool to find active 
black holes in this regime. Large samples of radio-loud active 
galatic nuclei (AGN) provide increasing evidence that two dif- 



* The master catalog (described in Table 1) and the catalog listing 
all radio matches per galaxy, are only available in electronic form at 
the CDS via http: //cdsweb.u-strasbg . fr/cgi-bin/qcat?J/A+ 
A/ or at http : //ragolu . science . ru . nl. 



ferent modes of accretion onto super-massive black holes indeed 
exist (Ghisellini & Celotti 2001; Kording et al. 2006; Ghisellini 
et al. 2011; Wu et al. 2011; Best & Heckman 2012). 

The first catalogs of extra-galactic radio sources were con- 
structed by gathering a (heterogeneous) set of optical follow- 
up observations (Schmidt 1968; Veron & Veron 1974; Kuehr 
et al. 1981; Laing et al. 1983; Jones & McAdam 1992). With 
the advent of wide-field optical surveys and deep radio surveys 
with high angular resolution, it has become possible to match 
cataloged sources at both wavelengths, and thus systematically 
construct large catalogs of radio sources with distance informa- 
tion (Condon et al. 2002; Ivezic et al. 2002; Best et al. 2005b; 
Sadler et al. 2007; Mauch & Sadler 2007; Kimball & Ivezic 
2008; Donoso et al. 2009; Brown et al. 2011). Such catalogs 
have played a key role in the study of star formation, AGN, and 
the interplay between them. 

While at intermediate to large redshifts (z > 0.1) one can 
survey a limited area of the sky to obtain a representative slice 
of the extra-galactic volume, studies of the local universe require 
full-sky coverage (i.e., 4n: solid angle) to map the anisotropy of 
the matter distribution. None of the modern catalogs of active 
black holes meet this requirement. 

With the recent release of the 2MASS redshift survey 
(2MRS; Huchra et al. 2012), currently the deepest all-sky red- 
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shift survey, it has become possible to construct an extra-galactic 
radio catalog that covers 90% of the volume of the local uni- 
verse. Indeed we made this our goal: obtain a complete, all-sky 
catalog of galaxies that emit at radio wavelengths. We have con- 
structed this catalog in the most systematic way possible to en- 
sure that the selections effects are well-defined (opposed to com- 
piling a list of known radio sources from the literature). Since the 
radio galaxies of our catalog are a subset of normal galaxies and, 
by construction, the sample is the largest of its kind, it present a 
powerful tool for a statistical study of the relation between black 
hole activity and galaxy environment. 

Our primary motivation for building an all- sky radio cata- 
log is to obtain a volume-limited sample of galaxies that could 
be powerful enough to accelerate ultra-high energy cosmic rays 
(UHECRs, charged particles with an energy in excess of 5 x 
10^^ eV) and study the magnetic field and energetics of these 
sources. For this study, full-sky coverage is key because inter- 
actions with photons of the cosmic microwave background limit 
the distance an UHECR can travel to about 100 Mpc, the so- 
called, GZK horizon (Greisen 1966; Zatsepin & Kuz'min 1966); 
if the sources of UHECRs are rare (< 10"^ Mpc~^), the full-sky 
has to be searched to obtain a sizable sample. 

A good candidate source of UHECRs is Cen A: this radio 
galaxy may just be powerful enough to accelerate protons up to 
the ultra-high energy scale (for a recent review see Biermann 
& de Souza 2012). Since Cen A is the nearest and one of the 
best studied radio galaxies, it provides a good anchor for the 
rest of the catalog; the flux of Cen A at 100 Mpc (~ 1 Jy at 
1 GHz and K = 11) is comfortably within the limits of existing 
surveys. Our main goal can thus be summarized in one sentence: 
find all radio galaxies within 100 Mpc that are as a luminous as 
Cen A. Since Cen A is relatively faint at optical wavelengths, this 
requirement actually allows us to find most typical radio galaxies 
within 200 Mpc. 

In this paper, we focus on the construction and properties of 
the radio-bright sample (a larger sample, obtained by lowering 
the radio flux limit, will be presented in a future publication). In 
section 2 we discuss our matching algorithm and we compare 
our sample to existing radio catalogs. In the third section we 
present the number counts (sec. 3.1), the cross-correlation of ra- 
dio galaxies with the local matter distribution (sec. 3.2), and the 
luminosity functions (sec. 3.3). We close with a discussion (sec. 
4). In the second paper in this series, we will discuss the mag- 
netic fields, jet power, and energy injection of the radio galaxies 
with particular emphasis of its relevance for UHECRs. 

2. Catalog construction 

As explained in the introduction, we wish to obtain a flux- 
limited, all-sky catalog of extra-galactic radio sources, from 
which a volume-limited sample can be derived. In this section 
we discuss the construction of this catalog. First, we will de- 
scribe the input (sec 2.1), followed by a detailed discussion of 
the matching algorithm (sec. 2.2). In section 2.3 we present the 
columns our catalog. We assess the completeness in section 2.4 
and probability for random matches in section 2.5. A compari- 
son to existing extra-galactic radio catalogs is presented in sec- 
tion 2.6. Finally, we briefly discuss some newly identified radio- 
emitting galaxies in section 2.7. 

2.1. Input 

In the following paragraphs we discuss the input for our cat- 
alog: the 2MASS Redshift Survey (2MRS), the NRAO VLA 
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Figure 1. The diff'erence between the Hubble distance, after our 
correction for peculiar velocity (/^Hubble), and the final distance 
(^finai)- The latter is the inverse- variance weighted mean of 
^Hubble and median redshift-independent distance from NED-D. 
The fractional uncertainty on the Hubble distance decreases with 
redshift, hence the convergence to zero diff'erence. The dashed 
line indicates the redshift cut that is applied in the computation 
of the luminosity functions and clustering (sec. 3). 

Sky Survey (NVSS), and the Sydney University Molonglo Sky 
Survey (SUMSS). 

2.1 .1 . Redshift survey 

To be able to select a volume-limited sample of radio sources, we 
need a measurement of the distance to each source. Currently, 
the deepest all-sky redshift survey is the 2MRS (Huchra et al. 
2012), which covers essentially the entire extra-galactic sky ac- 
cessible at optical wavelengths (91% of the full-sky). The tar- 
gets of this survey are galaxies from the 2 Micron All- Sky 
Survey (2MASS, Skrutskie et al. 2006), contained in the ex- 
tended source catalog (Jarrett et al. 2000; Jarrett 2004), selected 
with the following requirements: 

- < 1 1.75 and detected in the if -band 

- E{B-V)<\ 

- \b\ > 5° for 30° < / < 330°; \b\ > 8° otherwise. 

Where E{B-y) is the extinction from the maps of Schlegel et al. 
(1998); / and b are the Galactic longitude and latitude, respec- 
tively. Ks is the isophotal magnitude (measured in an elliptical 
aperture defined at the 20 mag per square arcsecond isophote). 
The ^-band is centered at /I ^ 2.2 yum, the near infrared (NIR) 
part of the electromagnetic spectrum. 

The 2MRS cuts select 44,599 galaxies from 2MASS; the 
current catalog contains redshifts for 97.6% of these. The me- 
dian redshift is 0.028, 90% of the sample is contained within 
z < 0.052. We applied no cuts on the 2MRS catalog. 

Some care has to be taken when converting redshifts of 
nearby-galaxies to distances using the Hubble law. We adopt the 
approach used by Blanton et al. (2005) for the construction of the 
NYU Value-Added Galaxy Catalog. First, we shift to the Local 
Group barycenter using the heliocentric velocity determination 
of Yahil et al. (1977). In this frame we estimate the most likely 
distance and its uncertainly using the model of the local veloc- 
ity field of Willick et al. (1997) based on the IRAS 1.2 Jy red- 
shift survey (Fisher et al. 1995). We adopt //q = 72kms"^ and 

= 0.3, Qa = 0.7 to convert the peculiar- velocity corrected 
redshifts to a distance. 
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For galaxies with a radial velocity (measured in Local Group 
barycenter) below 6, 000 km s"^ we attempt to improve the 
Hubble law distance using published redshift-independent dis- 
tances (e.g., Cepheids, Supernova Type la, Tully-Fisher rela- 
tion) as listed in NED-D^. For 3351 galaxies (8%) we com- 
pute the inverse variance weighted mean distance modulus from 
our peculiar- velocity corrected Hubble distance and the median 
redshift-independent distance modulus. The corrections to the 
luminosity distances derived solely from the Hubble law are 
modest, -0.0006 ± 0.07 dex in the mean. Corrections larger 
than 0.5 dex are observed only for galaxies with a Local Group 
barycenter radial velocity smaller than 187 kms"^ (Fig. 1). For 
the radio-emitting galaxies in our final sample (section 2.2.3), 
the coverage of NED-D is 100% below 250 kms"\ so the un- 
certainty on the radio luminosity due to peculiar velocities is not 
significant. 

2.1.2. Radio catalogs 

No single radio survey covers the full sky, we thus have to com- 
bine multiple surveys. To avoid sky-dependent selection eff'ects, 
these surveys have to be similar in depth, frequency and reso- 
lution. With these requirements only one set of radio surveys 
remains^: NVSS (Condon et al. 1998) and SUMSS (Bock et al. 
1999; Mauch et al. 2003), covering the Northern and Southern 
Hemisphere, respectively. 

The NVSS catalog is derived from images obtained with the 
VLA in D and DnC configurations at 1.4 GHz and covers the 
sky north of declination -40° (75% of the full-sky). The size of 
the point-spread-function (PSF) of these images is 45'' FWHM 
and the astrometric accuracy is better than 1" for > 15 mJy. 
The typical rms brightness fluctuations in the NVSS images are 
0.45 mJy/beam and the completeness limit is about 2.5 mJy. 

The SUMSS catalog is constructed from data obtained with 
MOST (Molonglo Observatory Synthesis Telescope; Mills 
1981; Robertson 1991) at 843 MHz. The resolution of the 
SUMSS images is 43" FWHM. The typical astrometric accu- 
racy in right ascension is 2 - 3" at 10 mJy and 1" above 20 mJy; 
position uncertainties in declination (S) are typically a factor of 
csc(l^l) higher than in right ascension. The rms noise level of 
the SUMSS images is ~ 1 mJy/beam, the completeness limit is 
about 10 mJy. The survey covers the sky south of declination 
-30° with \b\ > 10°. The sky below \b\ = 10° that is not in 
the NVSS footprint is covered by MGPS-2 (The second epoch 
Molonglo Galactic Plane Survey; Murphy et al. 2007) which is 
the Galactic counterpart to SUMSS. The MGPS-2 catalog, how- 
ever, currently only contains compact sources which makes it 
unsuited for our matching algorithm. We therefore removed the 
\b\ < 10 region for S < -40, reducing the area covered by our 
catalog from 91% to 88%. 

It is important to realize that both the SUMSS and NVSS 
catalog are constructed by fitting elliptical Gaussians to the radio 
maps. Each entry in the catalog is a separate Gaussian with a 
measured major/minor axis and integrated flux. Extended radio 
sources are often broken up into multiple Gaussians, hence a 
single astrophysical object is not necessarily represented by a 
single entry in the radio catalogs. 



^ http : //ned . ipac . caltech . edu/Library/Di stances/ 
(v5.1) compiled by L Steer and B. R Madore. 

2 Other radio surveys such as FIRST (Becker et al. 1995) or AT20G 
(Murphy et al. 2010), could be considered more suitable for matching to 
optical catalogs of galaxies. However this set of surveys does not meet 
our requirement of covering the entire sky at a similar frequency. 



2.1.3. Flux limit 

As explained in the introduction, our primarily goal to obtain 
a volume limited sample of radio galaxies as bright as Gen A. 
We use the flux of Cen A at the distance that contains 95% of 
the galaxies in 2MRS to find a conservative value for the radio 
flux limit. Using Fhoo = 1330 Jy for all the radio emission of 
Cen A (Cooper et al. 1965) at 1.4 GHz and a spectral index^ a = 
-0.6, we obtain the following flux limits: F^qq > 213, ^543 > 
289 mJy for NVSS, SUMSS. We apply these limits to the sum 
of the integrated radio flux of all Gaussians that are matched to a 
galaxy. Extending these flux limits to the completeness limits of 
the radio surveys (~ 10 mJy) will be the subject of future work. 

2.2. Matching 

The next step is to match the radio surveys to the redshift sur- 
vey to find radio-emitting galaxies. Our aim is to automate 
the matching as much as possible. This reduces (or at least 
parametrizes) human bias and substitutes expensive man-hours 
with cheap cpu-time. Since this is our first attempt at fully au- 
tomated cross-wavelength matching, we will also inspect all 
matches manually. 

Two challenges have to be met: (i) the FWHM of the PSF of 
the radio images is over an order of magnitude greater than the 
NIR images, (ii) many radio sources are resolved and will ap- 
pear in the radio catalogs as multiple entries which can be off'set 
from the galaxy by several arcminutes (e.g., giant radio galax- 
ies). Our approach to this problem is to proceed in two steps: 
first we match on the catalog-level (sec. 2.2.1), then we assess 
these matches at the image-level (sec. 2.2.2). We designed the 
first step to be "generous", i.e., all potential matches should be 
found, at the cost of a large background of false identifications. 
At the image-level this background is often trivially rejected. 

In a nutshell, our matching pipeline is based on the assump- 
tion that the galaxy is source or origin of the radio emission, but 
we allow this emission to be displaced and even disconnected 
from the galaxy; it is optimized to recover radio emission from 
both starforming galaxies, and systems with more complicated 
or asymmetric morphologies such as FR II galaxies or "head- 
tail" sources. 

2.2.1. Catalog-level matching 

We use the size of the radio source and the angular distance be- 
tween this source and the galaxy to define the first matching cri- 
terion: 

J,<MimXFWHM, OR di<dnm , (1) 

where di is angular distance between the ith radio entry and the 
galaxy center and FWHM/ is the deconvolved major axis of the 
radio Gaussian. We adopted the following cuts: = 3 and 
d\im = 2'. In section 2.4 we verify that these settings are indeed 
generous enough. 

For the coordinate matching we use kSmatch, a new^ and 
efficient implementation of 3 dimensional binary tree search. 
It can find matches between two sets of points on a sphere in 
0(Nlog(N)) time as opposed to the O(N^) time needed for a 
brute force search (R Schellart et al. 2012 in prep). 

^ We define the spectral index, a, by Fy oc y". 
^ kSmatch is available under GNU General Public License at 
pschella . github . com/k3match 
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NGC0383 UGC 07115 




RA (J2000) RA (J2000) 

Figure 2. Two examples to illustrate the challenges that arise when matching optical images with a resolution of ~ 2" to radio images 
with a resolution of 45". The contours show the radio emission from NVSS; 2MRS galaxies are labeled with green pentagons. The 
entries of the NVSS catalog are shown by purple squares, we see that both jets are decomposed into multiple elliptical Gaussians 
(6 for NGC 383 and 11 for UGC 7115). The radio entries matched to this galaxy by the angular distance requirement (Eq. 1) are 
labelled with crosses. For UGC 071 15, the remaining four jet components are included after applying a friends-of-friends algorithm 
(sec. 2.2.1). 



In some cases, the requirement set by Eq. 1 will miss a small 
fraction of the radio emission (e.g., UGC 7115 in Fig. 2: the 
start of the jet is matched to the galaxy, but the rest of the ra- 
dio emission, extending 6' further out, is not). This problem 
is solved by extending the matches using a friend-of-friend al- 
gorithm. We add a new radio entry to the group of matches if 
dij < max(A^iini x FWHM/, dn^), with dij the angular distance 
between the current radio match / and the potential new entry j 
(i.e., a link length equivalent to Eq. 1), and we repeat this proce- 
dure until no new matches are found. This approach successfully 
recovers the entire structure of nearly all extended radio sources 
(as demonstrated by the images in Appendix A). 

For each galaxy, the total radio emission is simply given by 
the sum of the integrated flux of all the components that are 
matched to the galaxy. Using Eq. 1 we obtain 1273 galaxies with 
a total radio flux above our flux limit (section 2.1.3), with a to- 
tal of 8452 matches between galaxies and entries in the radio 
catalogs. The generous criteria used at this stage are required to 
match galaxies with displaced radio counterparts, but they in- 
evitably yield false matches; we discuss the rejection of these 
false matches in the following section. 

2.2.2. Image-level rejection 

We now proceed to remove false matches from the set of 1273 
galaxies from the catalog-level matching (section 2.2.1) using 
the information contained in the pixels of the radio images. We 
stress that in nearly all cases these false matches are trivially 
identified (e.g., the left panel of Fig. 2, in this field all galaxies 
are matched to the radio emission, but the source of this emission 
is undoubtedly the central galaxy, NGC 0383). The only chal- 



lenge is to parametrize the judgement of a trained astronomer^. 
Again our approach is to use the galaxy as a starting point to 
identify (disconnected) radio emission as originating from a sin- 
gle source (e.g., IC 2722 in Fig. A.2). 

Our algorithm consists of the following steps (in Fig. 3 we 
show two examples). 

1. Make a cut-out of the radio image centered at the coordi- 
nates of the galaxy. The width of this image (/) is set by the 
maximum angular distance between the galaxy and the radio 
sources matched to this galaxy. 

2. Draw contours on this sub image. The level of the lowest 
contour is given by max((T(/)/7, 3 x (t(/5o)), where cr{I) is 
the standard deviation of the pixels of the image obtained in 
the previous step and /50 is subset of this image containing 
the pixels below the median of the image. 

3. Find all pixels contained by lowest contours that are within 
max(FWHM/, 30") of the location of the galaxy, where 
FWHM/ is the deconvolved major axis of the radio sources 
matched to this galaxy. If no pixels are found, we reject this 
galaxy. Most galaxies are rejected at this step since they are 
not connected to the radio contours (e.g., left panel of Fig. 
2). Radio sources that consist of two lobes (FR II morpholo- 
gies) are not rejected because the galaxy connects the two 
lobes (e.g., UGC 7360 in Fig. 3). 

4. Find all pixels within the lowest contour that are connected 
to pixels of the previous step. Any radio matches outside 
this group are rejected. This step removes unrelated radio 
emission that has been matched at the catalog-level (e.g.. 



^ This is a classical problem of Computer Vision; it may likely be 
solved more efficiently using tools from this field. 
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36-^ 30-^ 24' 18' n'n9"nT 36^ 30-' 24' 18' u'^imT 36' 30^ 24' 18' I2'n9'"12' 

RA (J2000) RA (J2000) RA (J2000) 




RA (J2000) RA (J2000) RA (J2000) 



Figure 3. The image-level matching algorithm at work (top: UGC 7360, bottom: NGC 0547). In step 1 & 2 we determine the size 
of the frame and draw contours. In step 3 we find the pixels above the value of the outer contour (labelled with red plusses) and 
in step 4 we link these pixels to the optical center of the galaxy (green squares). In step 5, we find the pixels that are connected to 
the group obtained in the previous step (shown by green circles), these define the radio emission matched to this galaxy. Finally we 
measure the geometrical center (white cross) and the flux- weighted center (cyan circled cross). For the latter, we use only the pixels 
that contain 80% of the total flux of the group (labeled with blue squares). 



the head-tail source that was initially matched to NGC 0547, 
Fig. 3). 

5. If more than one elliptical Gaussian from the radio catalogs 
is matched to a galaxy, we try to estimate the center of the 
radio emission. We compute the geometrical center by min- 
imizing 



(2) 



where dij is distance between pixel / and j and the sum runs 
over all pixels in the group (i.e., pixels obtained in step 4). 
We also compute a flux- weighted center by minimizing 



di = msix(dijWij) 



(3) 



where the weight Wfj is unity for the pixels that contain 80% 
of the total flux in the group and zero otherwise. The ge- 
ometrical method (Eq. 2) yields a robust measurement of 
the center for radio sources that have at least one symmetry 
axis that intersects the location of the galaxy, while the flux- 
weighted center gives a better estimate for radio sources with 
two bright spots at equal distance to the galaxy. 

Step 3 rejects 412 galaxies that are not close enough to the outer 
contour of the radio emission. After removing unrelated radio 
emission in Step 4, we al left with 769 galaxies above the radio 
flux limit. 



To define our final sample, we use a cut on the angular dis- 
tance of the galaxy to the geometrical and the flux-weighted cen- 
ter. Since these measures of the center are not alway appropriate 
(e.g., for head- tail sources), we also consider the minimum an- 
gular distance between the galaxy and the radio matches as well 
as the minimum angular distance between the galaxy and the 
second-highest contour (if the galaxy is within this contour we 
set this distance to zero). In Fig. 4 we show the distribution of 
these distances. Based on these distributions, we adopt the fol- 
lowing criteria for galaxies with multiple radio matches: 



d < 



1.2' to geometrical center 

1.2' to flux- weighted center 

15" to second highest contour 

8" to nearest radio catalog entry 



(4) 



Galaxies with multiple radio matches should pass at least one 
of these angular distance criteria. For galaxies with a single ra- 
dio match, the flux-weighted or geometrical center should not 
be used. Instead, we require at least one of the following two 
criteria: 



d < 



8" to radio catalog entry 
0" to second highest contour 



(5) 



The second of these two requirements implies that the coordi- 
nates of the galaxy are within the second highest radio contour. 
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Nearest radio Gaussian 
Single match 
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Figure 4. Histogram of angular distance between the NIR center of the galaxy (i.e., 2MRS coordinates) and different measure of 
the center for the radio emission for the candidate matches from the catalog cuts (Eq. 1). Left: the geometrical center (Eq. 2) and 
the flux-weighted center (Eq. 3). Right: the minimum angular distance between the galaxy and any of the matched radio sources for 
galaxies with a single match (thick blue line) and multiple matches (thin grey line). 



For point sources this requirement is superfluous (the galaxy al- 
ways lies within the highest contour). For extended radio emis- 
sion, however, the galaxy can be within the second highest con- 
tour while the distance to the radio catalog entry is > 8". 

The image-level cuts (Eqs. 4 & 5) reject 119 galaxies; 78 
sources are flagged for manual inspection because more than 
one galaxy is matched to the same radio source. Besides these 
78 ambiguous matches, we will also inspect the rejected and ac- 
cepted matches to measure how well our image-level algorithm 
is preforming. 

2.2.3. Manual inspection & classification 

The sample obtained after the catalog-level matching (Eq. 1) is 
small enough to inspect all matches. We stress that, at this point, 
we do not add new sources to the sample, the only goal is to iden- 
tify spurious or ambiguous results from the automatic matching 
pipeline. 

For 5 1 galaxies the radio emission matched to the galaxy by 
the image-level algorithm (sec 2.2.2) was rejected after human 
inspection, hence the purity of the sample after our fully auto- 
mated image-level matching procedure is 91%. Some of these 
initially accepted matches are very hard to classify or reject au- 
tomatically (e.g., blended sources or, very rarely. Galactic emis- 
sion), but most are due to artifacts in the radio images. We found 
that in only six cases, our image-level cuts (Eqs. 5 & 4) rejected 
a galaxy while at inspection we found this to be genuine match, 
implying an efficiency of 99%. 

The sources above the flux limit (sec 2.1.3) that remain after 
manual inspection comprise the final sample: 575 radio-emitting 
galaxies. Using the morphology of the radio emission, we clas- 
sify these into four categories, the number of sources in each 
class is given in brackets. 

- Point Sources (97): a single, unresolved radio source. 

- Starforming galaxies (52): extended radio emission that co- 
incides with the extended NIR emission from the galaxy. 

- Jets & Lobes (407): resolved radio emission beyond the NIR 
image of the galaxy that appears to originate from the cen- 
ter of galaxy. We did not attempt to subdivide this sample 
using the Fanaroff & Riley (1974) classification scheme be- 
cause most sources with an FR II morphology would not be 
classified as such based on their radio luminosity. 



- Unknown (19): if none of the above classes apply. These 
sources are not removed from our sample because they for- 
mally pass our image-level cuts (Eqs. 5 & 4). Yet their mor- 
phology suggests that the radio emission is not associated 
with the host galaxies. Users of the catalog may decide to 
remove the sources in this class. 

We stress that these classes are based on morphology only, they 
are not complete (i.e., some starforming galaxies or jets are clas- 
sified as point sources because the radio emission is not re- 
solved). 

2.3. Catalog columns 

We combine the properties of the 575 galaxies of our final sam- 
ple into a single master catalog; the columns of this catalog are 
described in Table 1. In Table 2 we list the 50 most luminous 
sources within 120 Mpc. A selection of the full catalog is shown 
in Appendix B. Since for some applications it is useful to have 
acces to the individual radio matches per galaxy, we also provide 
a basic catalog, listing all matched radio Gaussians per galaxy. 

2.4. Completeness 

The 575 galaxies of our master catalog are a subset of the poten- 
tial matches obtained from the catalog parameters (sec. 2.2.1). 
We can check if the latter sample was large enough by compar- 
ing the minium angular distance between the radio Gaussians 
and the galaxies to the cuts that we have adopted. In Fig. 5 we 
see that the catalog selection is not limiting the size of the fi- 
nal sample, because all galaxies are well within one of the two 
requirements placed by Eq. 1 . 

We also confirm that the small difference in frequency be- 
tween NVSS (1.4 GHz) and SUMSS (843 MHz) does not lead 
to a noticeable selection effect: the number of NVSS sources 
over the number of SUMSS sources is 3.65, which is very close 
to fraction of the area (within the 2MRS footprint) probed by 
these surveys, 3.74. 

Since we manually inspected all matches obtained from the 
catalog-level selection (sec. 2.2.1), the completeness of our mas- 
ter catalog is primarily determined by the completeness of the 
input catalogs. As discussed in section 2.1, the completeness 
of NVSS and SUMSS is nearing 100% at our flux limits and 
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Table 1. Columns of the master catalog 



Column name 


Units 


Description 


Source* 


2MASX 




Target name of the object from 2MRS database 


2MRS 


NED_id 




NED id of galaxy 


NED 


ra 


deg (J2000) 


Right Ascension of galaxy 


2MRS 


dec 


deg (J2000) 


Declination of galaxy 


2MRS 


1 


deg 


Galactic longitude 


2MRS 


b 


deg 


Galactic latitude 


2MRS 


Kmag 


mag (AB) 


^-band isophotal magnitude, corrected for Galactic extinction 


2MRS 


Kmag_err 


mag (AB) 


uncertainty on K 


2MRS 


z 




Heliocentric redshift 


2MRS 


zdist 




Peculiar- velocity corrected redshift 


2.1.1 


zdist_err 




Uncertainty on z_dist 


2.1.1 


D 


Mpc 


Inverse- variance weighted mean of distance from NED-D and peculiar velocity 
corrected Hubble distance (h = 0.72) 


2.1.1 


D_err 


Mpc 


Uncertainty on the former row 


2.1.1 


gal_type 




Galaxy morphological type code: -9 to 9, encodes Hubble sequence, 98 if galaxy 
has never been examined (Huchra et al. 2012, Table A8) 


2MRS 


n_nvss 


# 


Number of radio Gaussians matched to galaxy (after manual inspection) 


NVSS, sec. 2.2.3 


n_sumss 


# 


Number of radio Gaussians matched to galaxy (after manual inspection) 


SUMSS, sec. 2.2.3 


F1400 


mJy 


Sum of integrated flux at 1.4 GHz of all radio matches (zero if n_nvss=0) 


NVSS 


F1400_err 


mJy 


Uncertainty on the former row 


NVSS 


F843 


mJy 


Sum of integrated flux at 843 MHz of all radio matches (zeros if n_sumss=0) 


SUMSS 


F843_err 


mJy 


Uncertainty on the former row 


SUMSS 


Fsyn 


mJy 


Flux at 1.1 GHz, obtained from 1.4 GHz or 843 MHz flux using a = -0.6 


NVSS or SUMSS^ 


Lsyn 


erg s~^ 


vLy at 1.1 GHz 


NVSS or SUMSS^ 


sum_ma 


deg 


Sum of deconvolved FWHM of the major axis of all matched radio Gaussians 


NVSS or SUMSS^ 


inax_ma 


deg 


Largest deconvolved major axis FWHM of all matched radio sources 


NVSS or SUMSS^ 


lim_ma 


bool 


Limit flag on major axis (=0 if source is resolved) 


NVSS or SUMSS^ 


max_di st_to_gal 


deg 


Maximum angular distance between the galaxy and the radio matches 


NVSS or SUMSS 


inin_di st_to_gal 


deg 


Minimum distance between the galaxy and the radio matches 


NVSS or SUMSS 


contour_dist 


deg 


Angular distance to the second-highest radio contour (zero if galaxy is within 
contour) 


sec. 2.2.2 


geo_cen_ra 


deg 


RA of the geometrical center 


Eq. 2' 


geo_cen_dec 


deg 


Decl. of the geometrical center 


Eq. 2' 


flux_cen_ra 


deg 


RA of the flux- weighted center 


Eq. 3^ 


flux_cen_dec 


deg 


Decl. of the flux- weighted center 


Eq. y 


class 




Classification based on morphology: ^-Point Sources, q-Starforming galaxies, 
]-Jets & Lobes, \x-Unknown 


sec. 2.2.3 



Notes. * We lists what catalog provided the input for each parameter of our master catalog, we refer to a section or equation in this paper for 
parameters that not trivially obtained from the original catalog. In the overlap region (-40 < Decl. < -30) entries from both SUMSS and NVSS 
are matched to the same galaxy. How this information is combined depends on the parameter in question: ^ when both NVSS and SUMSS matches 
are available, we the pick the largest value. ^ when both NVSS and SUMSS data is available, we raise the major axis limit flag if the source is not 
resolved in both catalogs. ^ when both NVSS and SUMSS matches are available, we use the mean. 



the completeness of 2MRS is 97.6%. However, at the radio flux 
limit of our catalog, the systematic uncertainty on the total flux 
of sources with a complex radio morphology slightly decreases 
the completeness; a few sources may be missing from our final 
sample because our estimate of the total flux is below the true 
value of the total flux. This is not a serious limitation since most 
complex sources have been described in the literature which al- 
lows us to improve the estimate of their total flux (see section 
2.6). 

2.5. Probability of random matches 

To estimate the contamination due to background radio sources 
we shufiile the R. A. of the 2MRS galaxy catalog to obtain a ran- 
domized, uniform sample. The catalog-level cuts (Eq. 1) yield 
233 galaxies with a total of 1332 matches to the radio catalogs. 
Our image-level algorithm rejects 218 of these, leaving only 15 
galaxies. After checking this sample for imaging artifacts, we 
are left with 13 sources that genuinely passed the image-level 
cuts (Eqs. 4 & 5). Two of the 13 matches with the randomized 



2MRS catalog are radio sources that have already been unam- 
biguously identified as local radio-emitting galaxies using the 
original 2MRS data, hence these should not be included in the 
estimate of the number of background matches. 

All of the matches of the shufliled galaxy sample are to 
extended radio sources; most are background double sources 
whose geometrical center is within 1.2 arcminutes of the galaxy. 
As expected, they are classified into the Unknown morphological 
class (sec. 2.2.3). To conclude, the expected number of back- 
ground matches for extended radio sources in our catalog is 
^ 11, or 2% contamination. This number is of the same order 
as the number of galaxies in our catalog that have an unclassi- 
fiable radio morphology, confirming that most of these are very 
likely to be random matches. 

2.6. Comparison to other catalogs 

We now compare our sample of 575 radio-emitting galaxies with 
some well-known catalogues of extra-galactic radio sources. 
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Table 2. The 50 most luminous radio sources within 125 Mpc 



NED name vLy D Class 

(ergs~0 (Mpc) 




j 
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j 
j 
j 
j 
j 
j 
j 
j 
j 
j 
j 
j 
j 
j 
j 
j 
j 
j 
j 
j 
j 
j 
j 
j 
j 
j 
u 

j 
j 
P 
j 
j 
j 
P 
j 
j 
j 
j 
j 
P 
j 
j 
j 
j 
j 
j 
j 
P 



First, we consider the Kuehr et al. (1981) catalog, which is 
essentially complete for sources with F5000 > 1 Jy at 5 GHz and 
covers roughly the whole sky excluding \b\ < 10°. Three radio 
galaxies with z < 0.05 from this catalog are not found in our fi- 
nal sample: NGC 1265 (K = UA, z = 0.025), ESQ 252-GA018 
(K = 12.0, z = 0.034), and NGC 7385 (K = 13.3, z = 0.026). 
The ^-band flux of these galaxies is below the flux limit of 
2MRS (K < 11.75), hence by construction they are not included 
in our final sample. We shall not add these three galaxies to our 
final sample by hand because this would dramatically compli- 
cate the selection function — they would need to be excluded to 
obtain a well-defined volume-limited sample. We can also test 



Figures. The minium angular distance between the radio 
matches and center of the galaxy (di) and the minimum of this 
distance weighted by the major axis (J//FWHM/) for the 575 
galaxies of our master catalog. These parameters were used to 
select the sample of potential radio emitting galaxies (Eq. 1). We 
see that our cuts, di < max(2', 3 x FWHM/) (dashed red lines), 
are not limiting the number of galaxies in the final sample. 



the completeness of the Kuehr et al. (1981) catalog using our 
list of final sources. We limit our sample using the completeness 
limit of the Kuehr catalog F5000 > 1 Jy. converted to NVSS and 
SUMSS using a = -0.6 and restrict to at \b\ > 10, to find 16 
radio-emitting galaxies that are not listed in Kuehr et al. (1981). 

Next, we compare our catalog to that of Condon, Cotton, 
& Broderick (2002), obtained by a careful manual selection of 
UGC galaxies that are detected in NVSS. We restrict this cata- 
log using our ^-band and radio flux limits to 131 sources and 
find that five of these are not contained in our final sample. All 
five are nearby (z < 0.0017) starforming galaxies whose very ex- 
tended, low surface brightness radio emission has not been fully 
included in the NVSS catalog, causing them to fall (just) below 
our flux limit. For the main science goal of this work ("find all 
local radio galaxies"), these five galaxies are irrelevant. 

We also compare the total flux measured at 1.4 GHz by 
Condon et al. (2002) by a re-analysis of the NVSS images 
to our estimate using the sum of the fitted elliptical Gaussian 
listed in the NVSS catalogs. The flux they report is on aver- 
age 0.04 ±0.13 dex, higher than our work. Six radio galaxies 
show an off'set larger than 0.3 dex: NGC 0315, NGC 891, IC 342 
NGC 5127, NGC 6946, and NGC 7236. For NGC 0315, this is 
due a hotspot in the jet that has been missed by the friends-of- 
friends algorithm. The missing flux for the other sources is due 
to their complex radio morphology that is not entirely captured 
by the superposition of elliptical Gaussians. For these six cases, 
we use the total flux reported by Condon et al. (2002) as final 
value for the flux in our catalog, this reduces the the flux ofl'- 
set of the other 125 sources that are common to both catalogs to 
0.010 ± 0.06 dex. 

We should also compare our catalog to that of Jones & 
McAdam (1992), 193 southern extra-galactic radio sources with 
a flux density greater than 0.4 Jy and a size greater than 0.5'. 
After restricting this catalog to galaxies detected in 2MRS, we 
find three sources that are not in our final sample: WKK 4452, 
ESQ 137-G006 and ESQ 137-GOO. All three, however, are de- 
tected below \b\ < 10, which is outside the footprint of SUMSS 
(and thus outside the footprint of our catalog). In a future re- 
lease, we will complete our coverage of the southern sky below 



ESO 075- G 041 


41.80 


120 


NGC 1275 


41.37 


79 


2MASX J1713 1541-2502266 


41.26 


122 


UGC 01841 NOTESOl 


41.00 


98 


NGC 1128 NED02 


40.97 


97 


MESSIER 087 


40.92 


18 


NGC 5532 


40.92 


105 


NGC 3862 


40.90 


91 


Fornax A 


40.87 


21 


NGC 7236 


40.87 


116 


IC 4296 


40.84 


51 


2MASX J12543570-1234070 


40.79 


69 


CGCG 114-025 


40.77 


70 


PKS 0718-34 


40.77 


119 


UGC 11465 


40.76 


106 


NGC 0315 


40.73 


70 


NGC 6048 


40.63 


114 


NGC 0383 


40.63 


72 


NGC 6251 


40.60 


110 


ESO 460- G 004 


40.59 


105 


NGC 0547 


40.59 


77 


NGC 5090 


40.45 


47 


UGC 02755 


40.45 


108 


UGC 12064 


40.42 


74 


Cen A 


40.37 


3.6 


2MASX J17124278-2435477 


40.36 


103 


2MASX J01532586+71 15067 


40.35 


102 


NGC 1044 NEDOl 


40.31 


91 


2MASX J17204089-01 1 1573 


40.29 


124 


CGCG 403-019 


40.27 


119 


UGC 7360 


40.26 


32 


NGC 1167 


40.23 


74 


NGC 5127 


40.22 


71 


UGC 02783 


40.12 


91 


ESO 443- G 024 


40.11 


73 


NGC 6328 


40.11 


59 


UGC 11294 NOTESOl 


40.07 


67 


UGC 07115 


40.06 


96 


NGC 0741 


40.06 


80 


IC 0383 


40.04 


117 


2MASX J17354374-0720527 


40.02 


104 


NGC 4696 


40.00 


42 


NGC 3343 


40.00 


92 


NGC 6466 


39.97 


100 


NGC 7075 


39.96 


75 


UGC 12127 


39.96 


123 


NGC 2329 


39.95 


85 


2MASX J05444416+1648501 


39.95 


85 


NGC 5141 


39.95 


78 


IC 4374 


39.95 


90 
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\b\ = 10 using the MGPS-2 images and these radio galaxies will 
be included. 

We find a mean ratio between the 843 MHz flux measured 
by Jones & McAdam (1992) and this work of 0.015 ± 0.17 dex. 
Only Fornax A showed an ofl'set greater than 0.3 dex, which is 
due the complex morphology of this source. For this galaxy, we 
use the total flux measured by Jones & McAdam (1992) as final 
value for the flux in our catalog (this reduces the flux ofl'set to 
-0.014 ± 0. 1 dex). We also manually adjusted our catalog flux of 
Cen A to the total flux (i.e., inner and outer lobes) as measured 
by Cooper et al. (1965). 

Of the 575 radio sources in our final sample, 209 are not 
contained in the union of the three catalogs described above, 
plus the catalog of radio sources in the 6dFGS survey (Mauch 
& Sadler 2007) and the collection of known extra-galatic ra- 
dio sources and AGN of Veron & Veron (1974); Veron-Cetty 
& Veron (2006). 

2. 7. Newly found radio galaxies 

In the previous paragraph we found that over 30% of the galaxies 
in our sample are not contained in existing large-area samples of 
extra-galactic radio sources or AGN. Some of these have been 
identified as such by surveys that cover a smaller area of the sky, 
but in many cases this work is the first to classify these galaxies 
as radio emitting. Below we discuss some notable examples of 
these genuinely new identifications. We also list some sources 
that have been identified as (candidate) extra-galactic sources 
(e.g., based on their radio morphology or radio spectrum) by ear- 
lier work, but for which this work is the first to provide a redshift 
of the galaxy. The NED name of the galaxy can be used to find 
the images of these sources in Appendix A. 

- J00112171-h5231437 blended point source or head-tail 
source, member of a cluster that contains a radio relic (van 
Weeren et al. 2011). 

- NGC 0349 jet over 5' long, non-detection in the CRATES 
catalog of flat-spectrum radio sources (Healey et al. 2007). 

- J03204016-h2727485 complex geometry, multiple hotspots? 

- J03212595-hl806093 large (50 radio galaxy, shows both 
lobe emission and strong central emission. 

- NGC 1477 complex morphology, strong emission from 
galaxy center and faint lobes, contained in the CRATES cat- 
alog (Healey et al. 2007). 

- J05444416-h 1648501 extended jet emission, observed but 
not detected at VLBI scales: VERA 22 GHz flux < 0.1 Jy 
(Petrov et al. 2007). 

- J0612035 1-3257472 large source with a triple morphology 
(detected in both NVSS and SUMSS), galaxy is member of 
a rich cluster (Dressier 1980a). 

- J0733 1844-3654533 two large lobes (~ 10^ kpc), detected in 
both NVSS and MGPS-2. 

- Jl 1032753-4657471 large surface brightness difl'erence be- 
tween the two lobes. 

- IC 2722 "relaxed FR 11", looks like central engine has turned 
off'. 

- ESO 505- G 014 spectacular narrow and long (10') jet. 

- J13595556-1056266 clearly extended jet emission, detected 
by Chandra (Evans et al. 2010). 

- J16103572-0511173 jet morphology, but galaxy is not ex- 
actly at the center of the radio emission. 

- J 16390277-6505079 complex morphology, seems to show a 
hot spot displaced from current jet axis. Galaxy is located in 
the Great Attractor region (Radburn- Smith et al. 2006). 




Figure 6. Map in Galactic coordinates with the radio galaxies of 
the volume-limited sample (z < 0.03). The area of the circles is 
proportional to the radio flux of the source. The location of some 
famous sources is indicated (M87 and NGC 1275 are the bright- 
est members of the Virgo and Perseus cluster, respectively). 



- J17050125-2445099 two lobes and core emission, in 
Ophiuchus region (Hasegawa et al. 2000). 

- Jl 7 124278-2435477 complex morphology, in Ophiuchus re- 
gion (Hasegawa et al. 2000). 

- J17131541-2502266 two lobes, with yLy = 2x lO^^ ergs'^ 
this is one of brightest radio sources within 125 Mpc. It was 
identified as a potential double source (with unknown red- 
shift) by Slee, Siegman, & MulhaU (1982). 

- J19264081-h4123284 head-tail source. 

- J 19300 192- 1509191 resolved radio emission, most likely 
central jet emission. Detected by ROSAT (Voges et al. 1999), 
radio. X-ray and optical (without redshift) matched by Bauer 
et al. (2000). 

- J20430830-5059054 resolved radio emission, morphological 
class not clear. 

- J2 1492796-6429 194 slightly resolved, identified as radio 
source in the direction of the Abell cluster (i.e., without red- 
shift) by Robertson & Roach (1990). 

- J2238012 1-1-4 107363 two non-axisymmetric lobes. 

- MCG -05-55-032 wide-angle morphology, classified as a 
point source in the CRATES catalog (Healey et al. 2007). 

3. Analysis 

The 575 radio-emitting galaxies we identified in the previous 
section are a subset of a much larger flux-limited sample of 
"normal" galaxies, which allows for detailed comparison of both 
classes. In section we will work with two samples of radio galax- 
ies. We define sample A, the volume-limited sample that is used 
to measure the Hubble type abundance (sec. 3.1) and clustering 
(sec. 3.2) of radio galaxies. And we define sample B, which is 
limited only by the radio luminosity of the galaxy, to compute 
the luminosity functions (sec. 3.3). A detailed analysis of the 
magnetic field and jet power of the radio galaxies will be pre- 
sented in an accompanying paper (van Velzen et al. 2012). 

3.1. Number counts in a volume-limited sample 

Our master catalog is limited by both the flux limit of the red- 
shift survey and the radio flux limit we imposed. To obtain a 
volume-limited sample, we thus have to make a cut on both the 
^-band luminosity and the radio luminosity. We shall restrict 
our sample to z = 0.03, or comoving distance of 130 Mpc; the 
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Figure 7. The ^-band and 1 GHz luminosity (vLy) for all extra- 
galactic radio sources in our sample (top) and a volume-limited 
sample (bottom). Most sources that are morphologically classi- 
fied as AGN (i.e., the Jets and Lobes class) are more luminous 
than vLy = 10^^ erg s"^ Our volume-limited sample shows that 
the galaxy NIR luminosity and radio luminosity are not strongly 
correlated. 



corresponding minimum ^-band luminosity is Mk = -23.78, 
this limit is just high enough to select all galaxies with a K- 
band luminosity equal to Gen A (Mk = -23.8). Applying this 
distance and NIR luminosity cut to 2MRS yields 10964 galax- 
ies; our master catalog is reduced from 575 to 153 galaxies. The 
flux limit applied to the radio catalogs was F1400 > 213 mJy, 
^843 > 289 mJy for NVSS, SUMSS, respectively. Hence for 
z < 0.03, our catalog contains all sources with a radio luminos- 
ity greater than vLy = 5.8 x 10^^ erg s"^ at 1.4 GHz and vLy = 
4.8 X lO^^ergs-i at 843 MHz (or L1400 > 4.2 x 10^^ WHz'^ 
L843 > 5.6 X 10^^ WHz"^). Applying this radio luminosity re- 
quirement selects 461 galaxies of which 74 are below z = 0.03, 
all of these also obey the cut in ^-band luminosity. None of 
these radio-emitting galaxies are morphologically classified as 
Starforming. In the following we shall refer to this set of 74 
sources within z = 0.03 as sample A: the volume-limited sam- 
ple of radio galaxies. We list the properties of the galaxies in 
sample A in Table B.2. 

If we restrict sample A to galaxies with a radio luminosity 
that is greater or equal than Gen A (vLy = 2.6 x 10^^ erg s"^ at 
1 GHz), only 24 radio galaxies within z = 0.03 remain. 

In Fig. 6 we show a map in Galactic coordinates with the 
location of the 74 radio galaxies of sample A. The radio and K- 
band luminosity of both the full flux-limited sample and volume- 
limited sample are shown in Fig. 7. The 74 radio galaxies com- 
prise 0.65% of all galaxies in the volume-limited sample; in Fig. 



Radio galaxies (L1400 > 4 x 10^3 WHz"^) 
excluding Cen A 
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Figure 8. The fraction of radio galaxies in the volume-limited 
sample (sample A, defined in sec. 3.1) as a function of distance 
to the local group barycenter. The grey area shows the 90% CL, 
based on Poisson statistics for the number of observed radio 
galaxies. Since our sample probes over 90% of the volume, this 
uncertainty reflects the cosmic variance, rather than the statis- 
tical significance of this observation. Cen A remains the only 
radio galaxy up to 16 Mpc, where M 87 joins the ranks (at this 
distance, the number of normal galaxies in the volume-limited 
sample is 62). 



8 we show this fraction as a function of distance to the barycen- 
ter of the Local Group. 

Using the 2MRS morphological galaxy classification for 
sources with K < 11.25, we find the following Hubble type 
fractions for our volume-limited sample of radio galaxies: 47% 
Ellipticals, 47% SO galaxies, and 3% Spirals. For the ~ 10"^ nor- 
mal galaxies within z = 0.03 the E, SO and S fractions are 14%, 
31%, and 54%; for massive normal galaxies with Mk < -25, 
these fraction are 26%, 40%, and 33%. To allow a direct compar- 
ison to the galaxy morphology distribution of radio galaxies, we 
compute the Hubble type fractions of non-radio galaxies in bins 
of absolute ^-band magnitude and we compute the mean frac- 
tion weighted by the number of radio galaxies in each luminosity 
bin. We thus find the Hubble type fractions for non-active galax- 
ies drawn from the observed mass distribution of radio galaxies: 
27%, 38%, and 34% (E, SO, and S). 



3.2. Clustering of radio galaxies 

Our sample of radio galaxies is a subset of an all-sky redshift 
survey, allowing us to study the clustering of radio galaxies with 
respect to the matter distribution. We use the volume-limited 
sample of radio galaxies within z = 0.03 (sample A, defined 
in section 3.1) and we restrict to z > 0.003, yielding 73 

radio galaxies. In Fig. 9 we show a "top view" of this 
volume-limited sample; more than 50% of the radio galaxies are 
within the contour that encompasses half of the normal galaxy 
sample, an obvious sign of enhanced clustering. We quantify 
this clustering by counting the number of pairs between radio 
galaxies and normal 2MRS galaxies as a function of comoving 
distance, PRcid), i.e., the total galaxy count within a comoving 
distance d of each radio galaxy. We compare this to the num- 
ber of pairs for random subsets of 73 non-radio galaxies of the 
volume-limited sample, Pmatter(<^)- In the top panel of Fig. 10 we 
show Prg/ (Pmatter) - 1: a highly significant excess of clustering 
with respect to the normal galaxy distribution is clear. 
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Figure 9. Projected density of the volume-limited sample of ~ 
10"^ galaxies (z < 0.03) along /? = / = in Galactic coordinates. 
The contour encompass 95% and 50% of the sample. The lack 
for sources around j = is due to the disk of the Milky Way. One 
clearly sees that the radio galaxies in our volume-limited sample 
(blue circled crosses) reside in regions of high galaxy density. 

Clustering is well-known to correlate with galaxy mass and 
morphology (e.g., Dressier 1980b; Bamford et al. 2009) and the 
majority of the radio galaxies in our sample have massive, early- 
type hosts. Hence galaxy luminosity and morphology are con- 
founding factors in the analysis of radio galaxy clustering. One 
can remove the clustering due to host mass by computing the 
density excess with respect to the radio-loud galaxy mass dis- 
tribution. We thus measure (pmatter*) by drawing random sub- 
sets of non-radio galaxies from the observed ^-band luminosity 
distribution of radio galaxies. We show the result in the middle 
panel of Fig. 10; while the amplitude of the excess decreases, 
it remains highly significant. By forcing the random subsets of 
non-radio galaxies to have, on average, the same Hubble type 
(E, SO, or S) abundance as the radio galaxies we also correct the 
density excess for clustering due to morphology. We find that 
the galaxy number density around our sample of radio galaxies 
remains enhanced by a factor of 1.7 at 2 Mpc (Fig. 10 bottom 
panel); the probability that this enhancement is observed for a 
set of non-radio galaxies with the same luminosity and Hubble 
type distribution as the radio galaxies is less than 0.3%. 

The brightest galaxy of a cluster is often found to be radio- 
loud (Matthews, Morgan, & Schmidt 1964; Burns 1990) and this 
will contribute to the observed density excess of radio galaxies. 
We therefore preform a simple test to find the brightest cluster 
galaxies and estimate this contribution. First, we define a clus- 
ter or group as a collection of more than 5 galaxies from the 
volume-limited sample (Mk < 23.8) that can be connected by 
3 Mpc strings (i.e., a friends-of-friends cluster finding algorithm 
with a comoving link length of 3 Mpc). Of the 73 radio galaxies 
from sample A that we use here, 32 are in such groups. Of these, 
18 are the brightest member of their group (or cluster). After 
removing these galaxies, the enhanced clustering (corrected for 
mass and morphology of the host galaxies) of the remaining 55 
radio galaxies falls between the 2(T-level and the 3cr-level for all 




10° 10^ 

Distance (Mpc) 



Figure 10. The fractional density excess of radio galaxies in the 
volume-limited sample (defined in sec. 3.1). Top: the excess 
with respect to a volume-limited sample of non-radio galaxies 
(i.e., z < 0.03 and Mk < -23.78). Middle: the excess with 
respect to normal galaxies drawn from the ^-band luminosity 
function of radio galaxies. Bottom: the fractional density excess 
corrected for clustering induced by the morphological type of the 
host galaxies and their luminosity. Even after correcting for the 
clustering of their massive and predominantly early-type hosts, 
radio-loud jets are observed to occur in regins of high galaxy 
density. 



distances we consider; at 2 Mpc, the fractional density excess is 
50%. 

3.3. Luminosity function 

In this section, we use the 1/Vmax method (e.g., Schmidt 1968), 
to obtain the number density of radio galaxies and normal galax- 
ies as a function of luminosity. The l/Vmax method allows one 
to compute this density without applying a hard distance limit 
(such as the z < 0.03 limit that was used to define sample A in 
section 3.1). In this section, we define sample B: a set of power- 
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Galaxy luminosity {Mk — 5 log^g ^) Galaxy luminosity (log^o^/^o) 

Figure 11. Left: the ^-band luminosity function for powerful radio galaxies (sample B, L > 10^"^ WHz"^) and normal galaxies. We 
also show Schechter function fits to galaxies from the UKIDSS Large Area Survey (Smith et al. 2009) and the 6dF Galaxy Survey 
(Jones et al. 2006). The number density of powerful radio galaxies peaks at an absolute ^-band magnitude of -25.5 (using h = 0.72). 
Right: the fraction of radio galaxies as a function of galaxy luminosity. We see at rapid increase from zero for L < 5 x 10^^ L© to 
~ 20% for the brightest galaxies. We also show the fraction of radio-loud AGN (with the same minimum radio luminosity) in SDSS 
(Best et al. 2005a), converted to the ^-band luminosity using a mean mass-to-light ratio M/Lk = 0.8 (Bell et al. 2003). 
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Figure 12. The radio luminosity (vLy) function of extra-galactic radio sources at diff'erent frequencies. Right: 1.4 GHz from NVSS 
(S > -40) and 843 MHz from SUMSS (5 < -30, \b\ > 10). Left: the 1.1 GHz sample includes sources from both surveys by scaling 
to this frequency using the mean spectral index of a = -0.6. For vLy > 10^^ ergs"\ the luminosity function for radio galaxies 
with a jet or lobe morphology (open diamonds) is similar to full sample, while below this scale the contribution from starforming 
galaxies becomes clear. 



ful radio galaxies with L1400 > 10^^ WHz" . As before, we also 
require a minimum redshift z > 0.003 and exclude the sources 
with unknown morphology, because these are most likely ran- 
dom matches. The luminosity and minimum redshift cut have 
been picked to aid a comparison to previous work. These cuts 
leave 527 radio-emitting galaxies, with a median redshift of 
0.042, and 401 powerful radio galaxies with a median redshift 
of 0.052. 

We find (V/Vmax) = 0.494 ± 0.012 and (y/Vmax) = 0.496 ± 
0.014 for all radio-emitting galaxies and the powerful radio 
galaxies, respectively, which shows that density evolution with 
redshift or distance-dependent selection eff'ects are not a signifi- 
cant influence on the derived radio source densities. 

In Fig. 11 we show the ^-band luminosity function. One 
clearly sees that the powerful radio galaxies reside in the most 
luminous hosts and their density peaks at Mk ~ -25.5. The 
fraction of radio galaxies is a strong function of galaxy lumi- 
nosity. No radio galaxies with L1400 > lO^'^WHz"^ are found 
for Mk > -23.7 and the fraction rises rapidly to an asymptotic 
value of 0.25 for the most massive galaxies. We also compute 



the luminosity function of normal galaxies which agrees rea- 
sonably well with Schechter functions that have been derived 
form the UKIDSS Large Area Survey (Smith et al. 2009) and 
the 6dF Galaxy Survey (Jones et al. 2006) (except at the bright 
end, which is a known issue with Schechter functions). 

Recalling that the radio-emitting galaxies in our catalog orig- 
inate from NVSS (at 1.4 GHz) and/or SUMSS (at 843 MHz), 
we first compute the luminosity function for each survey sep- 
arately (Fig. 12, left panel). We find that, as expected, the lu- 
minosity functions from NVSS and SUMSS are nearly identi- 
cal. To increase the statistical power, we combine both surveys 
by scaling the flux to 1.1 GHz using the mean spectral index 
of a = -0.6 (Fig. 12, right panel). We find that the luminosity 
function of radio-emitting galaxies with a jet or lobe morphol- 
ogy agrees well with the results of the entire sample, except be- 
low vLy = 10^^ erg s"^. Below this luminosity, the contribution 
of starforming galaxies starts to dominate (e.g., Condon et al. 
2002). At 10^^ ergs"^ our estimate of the number density is a 
factor 2-3 above the 1 .4 GHz luminosity function derived form 
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2661 radio-loud AGN at 0.003 < z < 0.3 (Mauch & Sadler 
2007). Below, we suggest an explanation for this difference. 

4. Conclusions & Discussion 

Our conclusions can be summarized as follows. 

1. A new extra-galactic radio catalog covering 88% of the ce- 
lestial sphere has been constructed. The catalog contains all 
galaxies brighter than K = 1 1.75 that show radio emission at 
^1400 > 213 mJy or ^843 > 289 mJy; it contains 575 sources. 
Over 30% of these radio-emitting galaxies are not contained 
in existing large-area catalogs. 

2. Our matching algorithm identifies extra-galactic radio 
sources with 99% efficiency and 91% purity. All matches 
have been inspected manually to remove imaging artifacts 
and blended sources. The contamination due to background 
radio sources in the final sample is 2%; the Unknown mor- 
phological class contains 19 potential random matches. 

3. For a volume-limited sample within z = 0.03 (130 Mpc), we 
find 27 galaxies with a radio and ^-band luminosity greater 
or equal to Cen A, which is 0.2% of all 2MRS galaxies 
within the same volume. 

4. No obvious correlation between galaxy luminosity and total 
radio power is observed (Fig. 7). 

5. The ^-band luminosity function of radio galaxies peaks at 
Mk = -25.5 or a luminous mass of ~ 10^^ M© (Fig. 11). 

6. Of radio galaxies within z = 0.03, 94% are of the E/SO 
Hubble type; this higher than the fraction of bright (Mk < 
-25) normal galaxies that are of this morphological class 
(66%) or the E/SO fraction of non-active galaxies scaled to 
the same /sT-band luminosity distribution as the radio galax- 
ies (65%). 

7. The fraction of radio galaxies as a function of ^-band lumi- 
nosity rises asymptotically to 20% (Fig. 11). 

8. We have computed the fractional density excess of radio 
galaxies in a volume-limited sample as a function of dis- 
tance from the radio galaxies. Radio galaxies are signifi- 
cantly more clustered than normal galaxies; in a sphere of 
2 Mpc centered on the radio galaxies, the galaxy density is 
2.7 times greater than around normal galaxies (Fig. 10 top 
panel). 

9. After correcting the latter result for the extra clustering in- 
duced by the mass and Hubble type of the radio galaxies, 
the density excess remains significant (> 3cr): at 2 Mpc the 
density around radio galaxies is 1.7 times higher than around 
non-radio galaxies with the same mass and galaxy morphol- 
ogy distribution (Fig. 10 bottom panel). 

The observation that radio-loud AGN are found in the most 
massive galaxies is well-known. Indeed, the median absolute K- 
band magnitude of radio galaxies in the NVSS-6dFGS sample 
(Mauch & Sadler 2007) of Mk = -25 A is in excellent agree- 
ment with our results. The radio luminosity where the density of 
AGN overtakes that of starforming galaxies (vLy = 10^^ erg s"^) 
agrees with radio luminosity functions from other surveys at 
similar redshift (Condon et al. 2002; Mauch & Sadler 2007). 

The decreasing fraction of radio galaxies with distance (Fig. 
8) appears at odds with the cosmological trend of increasing 
AGN activity to z ~ 2 (e.g., Hopkins et al. 2007), yet most 
observers in other galaxies would find the same. Because radio 
galaxies are found in high density regions and the average as- 
tronomer lives in a galaxy inside a cluster, the local fraction of 
radio galaxies is typically observed to be higher than the cosmic 
average. 



For a standard mass-to-light ratio in the ^-band, M/Lk ~ 1, 
we find that the fraction of galaxies with a radio luminosity 
greater than L= 10^^ W Hz"^ is a factor 5 larger than the fraction 
of AGN with the same minimum radio luminosity in SDSS (Best 
et al. 2005a, 2007). We also find a factor 2-3 higher number 
density of radio sources at vLy = 10^^ ergs"^ (or 10^^ WHz"^) 
compared to the NVSS-6dFGS sample (Mauch & Sadler 2007). 
If we restrict our sample to sources detected beyond 200 Mpc 
(leaving 237 radio galaxies), while keeping V^ax of each source 
fixed, these two discrepancies disappear. Hence a large-scale lo- 
cal overdensity of radio galaxies or a distance dependent selec- 
tion eff'ect are the most likely explanations for the larger radio 
galaxy density derived from our sample. We note that Mauch 
& Sadler (2007) find (V/Vmax) = 0.532 ± 0.006 for their radio- 
loud AGN; this significant off'set from 0.5 implies evolution (or 
selection bias) is present in their sample. 

The conclusion of Ledlow & Owen (1995) that the number 
of radio galaxies detected in Abell clusters simply scales with 
the number of galaxies surveyed, appears incompatible with our 
detection of enhanced clustering of radio galaxies with respect to 
normal galaxies. Our clustering detection, however, is in agree- 
ment with more recent work at intermediate redshift (z > 0.1). 
Using the SDSS galaxy sample, Kauff'mann et al. (2008) found a 
factor ^ 2 increase for the galaxy counts in a projected radius of 
300 kpc around radio-loud AGN compared to radio-quiet AGN; 
Donoso et al. (2010) found that the projected cross-correlation 
at 2 Mpc between radio-loud AGN and luminous red galaxies 
(LRG) is 25% stronger than the auto-correlation of the radio- 
quiet LRG. Tasse et al. (2008a) compared the environment of the 
host galaxies of powerful radio sources (L1400 ^ 10^^ WHz"^) 
in the XMM-LSS field (Tasse et al. 2008b) to normal galax- 
ies of the same mass and redshift. For M > 10^^-^ M© their 
photometric-redshift-based density estimator finds an overden- 
sity at 450 kpc, close to the distance where the enhanced cluster- 
ing of radio galaxies in our sample exceeds the 3(T-level. 

The enhanced density around radio galaxies suggests a 
causal relation between galaxy environment and jet power. Since 
this excess is measured with respect to galaxies of the same mass 
and morphological type, the mechanism behind this relation is 
not primarily driven by black hole mass. The gas cooling out of 
the hot atmospheres of the host galaxies is often suggested as a 
potential mechanism to turn on radio-loud jets (Burns 1990; Best 
et al. 2005a). An interesting future application is to compare the 
luminosity of the core of the jet, which will be measured by a 
VLBA 8 GHz survey of 2MASS galaxies (Condon et al. 2011), 
to the large scale radio emission contained in our catalog. This 
will allow a study of jet power as a function of black hole mass 
— the latter can be estimated using the luminosity and Hubble 
type of the galaxy (e.g., Caramete & Biermann 201 1). 

The automated image-level matching presented here could 
be considered an improvement to the manual classification that 
has been used to construct similar samples in the past. For fu- 
ture studies of existing extra-galactic catalogs (e.g., the 2MASS 
Extended Source Catalog) as well as upcoming radio surveys 
(e.g., LOFAR, ASKAP), automated cross- wavelength identifi- 
cation will be key to manage the large number of sources. 
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Appendix A: Example images 

Here we show images of examples of newly identified radio- 
emitting galaxies discussed in section 2.7. For the images 
of all 575 radio-emitting galaxies in our sample we refer to 
the online catalog, http://ragolu.science.ru.nl/hcat. 
html. The following legend applies to all images: 

- Green pentagram: 2MRS galaxy. 

- Green cross (x): center of the frame (i.e., galaxy in question). 

- Purple square (□): SUMSS/NVSS catalog source. 

- Purple cross (x): matched catalog source by image-level al- 
gorithm. 

- Purple plus (-h): manually accepted match. 

- Large magenta cross: geometrical center of radio emission. 

- Large red circled cross: flux- weighted center of radio emis- 
sion. 



14 



Sjoert van Velzen et al.: Radio galaxies of the local universe 



2MASXJ001 12171+5231437 (DSS2-IR) 




2MASX JOOl 12171+5231437 (NVSS) 



NGC 0349 (DSS2-IR) 



NGC 0349 (NVSS) 



36' 28' 



RA (J2000) 

2MASXJ032040 16+2727485 (DSS2-IR) 



RA (J2000) 

2MASXJ032040 16+2727485 (NVSS) 




'0' 



RA (J2000) 

NGC 1477 (DSS2-IR) 




RA (J2000) 

2MASX J0612035 1-3257472 (DSS2-IR) 





RA (J2000) 

2MASX J0612035 1-3257472 (NVSS)_ 




RA (J2000) 

2MASX J0733 1844-3654533 (DSS2-red) 



RA (J2000) 

2MASXJ0733 1844-3654533 (MGPS-2) 






24-' 02"'08' 52' 

RA (J2000) 



2MASX J032 12595+ 1806093 (DSS2-IR) 



RA (J2000) 

2MASX J032 12595+ 1806093 (NVSS) 




2MASXJ054444 16+ 1648501 (DSS2-IR) 




RA (J2000) 

2MASXJ054444 16+ 1648501 (NVSS) 



48' 40' 5'"44'"32' 

RA (J2000) 



2MASX J0733 1844-3654533 (DSS2-red) 




RA (J2000) 

2MASX J0733 1844-3654533 (NVSS) 






2MASXJ1 1032753-4657471 (DSS2-IR) 



RA (J2000) 

2MASXJ1 1032753-4657471 (SUMSS) 



20' 33'"00' 7'<12'"4 

RA (J2000) 



34"'oo' 40' 20' ^^rm i^^i^vf 
RA (J2000) 




20' ii'i()3" 

RA (J2000) 



Figure A.l. Examples of newly identified radio-emitting galaxies. 



15 



Sjoert van Velzen et al.: Radio galaxies of the local universe 



IC 2722 NEDOl (DSS2-red) 



IC 2722 NEDOl (NVSS) 




RA (J2000) RA (J2000) 

2MASX J13595556-1056266 (DSS2-red) 2MASX J13595556-1056266 (NVSS) 




4'>00'"00' 56' 52-' 48' 13''59"'4. 

RA (J2000) 



, 2MASX 16390277-6505079 (DSS2-IR) 



RA (J2000) 

2MASX 16390277-6505079 (SUMSS) 




RA (J2000) 

2MASX 17124278-2435477 (DSS2-red) 





-24°30' 


^ 1 \ \ \ 

o 




33' 








1 










42' 




13-00' 48' 36' 17'42'"24' 1 

RA (.I2()()()) 


n"'m 48' 36' 17''12"'24' 

RA (J2000) 




ESQ 505- G 014 (DSS2-IR) 




28' 07"' 12' 56' 40' 

RA (J2000) 



ESQ 505- G 014 (NVSS) 







28' 07" 12' 56' 40' 12*06" 

RA (J2000) 



2MASX J16103572-0511173 (DSS2-IR) 2MASX J16103572-0511173 (NVSS) 

-5°08' ! 




RA (J2000) 

2MASX 17050125-2445099 (DSS2-IR) 




RA (J2000) 

2MASX 17050125-2445099 (NVSS) 




RA (J2000) 

2MASX 17131541-2502266 (DSS2-red) 




RA(J2000) 

2MASX J19300192-1509191 (DSS2-IR) 
,. i-"'-.Lf.. - / ' . ' • 

-15°07'h < 




RA (J2000) 

2MASX 17131541-2502266 (NVSS) 



RA (J2000) 

2MASX J19300192-1509191 (NVSS) 



04' 30"'00' 56' 19''29""5: 

RA (J2000) 



04' 30"'00' 56' 19''29'"52' 

RA (J2000) 



Figure A.2. Examples of newly identified radio-emitting galaxies. 
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Figure A.3. Examples of newly identified radio-emitting galaxies. 
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Appendix B: Catalog for the volume-limited sample 

Below we present a selection of the columns of the master cata- 
log for the radio galaxies of volume-limited sample (sample A, 
defined in sec. 3. 1). The entire catalog is available in CDS. It can 
also be obtained at http : //ragolu . science . ru . nl. 
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Table B.2. Volume-limited sample of radio galaxies. 
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Table B.2. continued. 
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fJv) 


fJv) 


Qo2iA er2 s~^) 


(arcmin) 




20520232-5704076 


313.00974 


-57.06877 


0.01175 


46.6 


8.8 




2.2 


39.72 


1.7 


p 


21313299-3837046 


322.88742 


-38.61790 


0.01817 


75.0 


9.7 


0.98 


1.1 


39.92 


2.6 


j 


21570595-6941236 


329.27487 


-69.68994 


0.02827 


119.8 


10.3 




44 


41.85 


5.6 


22041766+0440021 


331.07358 


4.66722 


0.02700 


119.4 


11.2 


0.78 




40.23 


2.7 


j 


22144500+1350476 


333.68747 


13.84653 


0.02616 


116.2 


10.3 


3.3 




40.83 


6.1 


22312062+3921298 


337.83585 


39.35825 


0.01713 


73.7 


9.9 


2.9 




40.38 
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22382946+3519400 


339.62277 


35.32798 


0.02759 


123.4 


9.7 


0.36 




39.92 


0.52 




22492263-3728185 


342.34439 


-37.47178 


0.02867 


124.0 


10.1 


0.34 


0.47 


39.91 


0.97 


j 
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